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Abstract. The quest for atmospheric spectral signatures that may witness biological ac- 
tivity in exoplanets is focused on rocky planets. The best targets for future, challenging 
spectroscopic observations will be selected among potentially habitable planets. Surface 
habitability can be quantified and explored with climate and atmospheric models according 
to temperature-based criteria. The conceptual, modellistic, technological and interpretative 
complexity of the problem requires to develop flexible climate and atmospheric models 
suited for a comprehensive exploration of observationally unconstrained parameters, and 
to simulate and interpret definitely non-terrestrial conditions. We present a summary and 
preliminary results on the work we are performing on multi-parametric explorations of the 
habitability and observational properties of rocky planets. 
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1. Introduction 


The surface environmental conditions that may 
allow a planet to be tagged as potentially hab- 
itable is the outcome of a plethora of physi- 
cal factors involving all scales — at the plane- 
tary, orbital, and stellar level, and continuously 
acting since the same process of planet for- 
mation. Therefore both the search for poten- 
tial habitable candidates as best spectroscopic 
targets, as well as the ensuing challenging 
task of interpreting observed data, possibly in 
terms of biosignatures, rest on thorough theo- 


retical multi-parametric physical modelling of 
the planetary surfaces and atmospheres. A vast 
range of unknowns affects also our understand- 
ing of the early evolution of Earth, when proba- 
bly terrestrial life was able to originate, as well 
as of Mars and Venus, that underwent com- 
pletely different evolutionary paths. We have 
developed climate and atmospheric modelling 
tools specifically tailored to perform habitabil- 
ity and interpretative studies, to be also ap- 
plied to simulate future space missions. Here 
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we present a short summary with preliminary 
results of our work on these subjects. 


2. Climate and atmospheric model 


In order to cast light on the physical properties 
and habitability of rocky exoplanets, as well 
as on the early evolutionary phases of Earth, 
Mars and Venus, we need to run climate sim- 
ulations. The choice of the models depends 
on how many observational data are available 
for individual cases. In general, for exoplan- 
ets we can measure only a few structural (ra- 
dius, mass) and orbital (semi-major axis, ec- 
centricity) parameters, plus a few properties of 
the host star (luminosity, spectral type, chem- 
ical composition). Spectroscopic data of the 
thin atmospheres of rocky exoplanets are still 
quite scarce and are limited to planets in tight 
orbits around very cold stars, conditions that 
are not particularly favourable for habitabil- 
ity. Given the lack of physical/chemical ob- 
servational constrains on key quantities, cli- 
mate/atmospheric models used in this field 
must be able to physically handle them as 
free parameters, in contrast with models used 
for the modern Earth, where most parame- 
ters are embedded. Here we briefly describe 
a model of this type that has been developed 
at INAF-Trieste in collaboration with clima- 
tologists based at CNR-IGG Pisa and Torino 
Politecnico. 


2.1. The climate model: ESTM 


The Earth-like planet Surface Temperature 
Model (ESTM) was developed to quan- 
tify the surface planetary habitability start- 
ing from the spatial-temporal distribution of 
surface temperatures (Vladilo et al. 2015). 
The core of ESTM is an energy balance 
model (EBM) with time-dependent quantities, 
complemented by: radiative-convective atmo- 
spheric column calculations, a set of phys- 
ically based parametrizations of the merid- 
ional transport, and descriptions of surface and 
cloud properties more refined than in clas- 
sic EBMs. Due to its dependence on time, 
latitude and vertical atmospheric column, the 
ESTM is in practice a seasonal, 2D climate 


model. Thanks to its extremely low computa- 
tional cost, it can be used to explore the param- 
eter space unconstrained by observations of in- 
dividual rocky (exo)planets, or to perform sta- 
tistical studies of climate properties whenever 
many runs for different parameter configura- 
tions are needed. The model has been recently 
upgraded and its most recent version, EOS- 
ESTM (Biasiotti et al. 2022), includes a new 
code, EOS, (Simonetti et al. 2022; Simonetti 
2022) for the treatment of the radiative trans- 
fer, as well as a number of parametrizations 
to simulate the climate impact of oceans, land, 
ice, and clouds as a function of temperature 
and host star zenith distance. EOS-ESTM can 
be applied to a large variety of rocky planets, 
with terrestrial and non-terrestrial atmospheric 
compositions illuminated by solar- and non- 
solar-type stars. One limitation of this (and any 
EBM- based) model is that it cannot be applied 
to tidally-locked planets. The model is under- 
going a continuous process of cross-validation 
with a variety of climate models with different 
levels of complexity. 


2.2. The atmospheric model: EOS 


EOS (named after the Greek goddess of Dawn) 
is a radiative-convective equilibrium model 
for single atmospheric columns. It is derived 
from the opacity calculation tool HELIOS-K 
(Grimm & Heng 2015; Grimm et al. 2021) and 
the radiative transfer code HELIOS (Malik et 
al. 2017; 2019), of the Exoclime Simulation 
Platforn]!] and has been specifically developed 
to model the relatively cold and dense at- 
mospheres of terrestrial planets. Special care 
has been devoted to collision-induced absorp- 
tion (CIA) continua and to the impact of the 
moist convection on the tropospheric verti- 
cal pressure-temperature profile. As its parent 
codes, EOS is GPU-accelerated and thus com- 
putationally efficient. 

EOS operates in the two-stream approxi- 
mation (see e.g. Pierrehumbert 2010) includ- 
ing non-isotropic scattering. Opacity calcula- 
tions can be carried on either line-by-line or 
by using the k-distribution method (Goody & 


! https://github.com/exoclime 
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Yung 1989). The code can take as input any 
line list or CIA table from popular spectroscop- 
ical repositories such as HITRAN (Gordon et 
al. 2022), HITEMP (Rothman et al. 2010) and 
ExoMol (Tennyson et al. 2020) and does not 
depend on hard-wired opacity parameteriza- 
tions. Similarly, the incoming stellar radiation 
spectrum can be specified by the user or set 
to an ideal blackbody of any temperature. The 
ability to change freely both the atmospheric 
composition and the instellation makes EOS 
extremely flexible and easy to keep updated. 

The primary output of EOS is the full re- 
flectance/emission spectrum for a specified set 
of parameters. From this we get the integrated 
quantities required by the EBM equation (the 
Outgoing Longwave Radiation -OLR-, and the 
top-of-atmosphere -TOA- albedo). The syn- 
thetic spectra allows to connect the climate 
state and habitability, as calculated by ESTM, 
to observations that will be performed with fu- 
ture instruments (Sec. [4). 

In practice, radiative transfer calculations 
for our climate simulations are performed by 
following the reverse procedure of Kasting (see 
e.g. Kasting et al. 1993). For any given atmo- 
spheric composition we obtain the OLR and 
the TOA albedo as a function of the surface 
temperature and reflectivity, the stellar zenith 
distance and potentially other variables, stor- 
ing them in lookup tables. These tables can in 
principle be used also by other climate models, 
thus being another standalone product of EOS- 
ESTM. 


3. Model applications 
3.1. Identifying habitable Exoplanets 


Most of the currently known exoplanets have 
been uncovered with indirect detection tech- 
niques, i.e. the transit and the radial velocity 
methods, in which the planet's influence on 
the host star is observed — luminosity dim- 
ming during transits, or Doppler shift of spec- 
tral lines due to movement by reciprocal grav- 
itational interaction. Therefore, larger planets 
in close orbits around small late-type stars, able 
to determine stronger planet-to-star signals, are 
preferentially selected (e.g. Winn 2018). Due 


to these observational biases, currently only 
one potential terrestrial analogue orbiting the 
habitable zone (HZ) of a G-type star has been 
detected, Kepler-452b (Jenkins et al. 2015). 

In Silva et al. (2017a) we constrained its 
current and past habitability for different atmo- 
spheric pressures and compositions, also ac- 
counting for the luminosity evolution of the 
host star. We found that a relatively low CO, 
abundance (x 0.04 bar) would allow a cur- 
rent habitable surface. Instead, a long-term > 2 
Gyr habitability, that could promote formation 
and widespread persistence of biological activ- 
ity able to pollute the atmosphere, would have 
required 2 100 times more greenhouse effect at 
early times due to stellar evolution effects (as 
envisaged also for the habitability of the paleo- 
Earth, Sec.3.2). 

The observational (and statistical) bias 
against detecting potentially habitable rocky 
planets is partly tackled for exoplanets or- 
biting the HZ of M-type stars. Even though 
the structural planetary properties could be of 
terrestrial-type (as e.g. the Trappist planetary 
system), it is expected that the star-planet in- 
teraction and evolution determines radically 
non-terrestrial situations, among which is the 
possible establishment of the spin-orbit tidal 
lock. This configuration is not adequate for 
EBM-type models (but see e.g. Haqq-Misra 
& Hayworth 2022). On the other hand, suf- 
ficiently eccentric orbits can severely delay 
or avoid the onset of the synchronous spin- 
orbit rotation. One such case is Gl-514b, a 
super-Earth recently detected by Damasso et 
al. (2022). We are investigating with EOS- 
ESTM the habitability of this planet as a func- 
tion of a large range of atmospheric properties. 
Very preliminary calculations show that indeed 
habitable solutions do exist, with possibly large 
seasonal excursions (Biasiotti et al., in prep.), 
that will consequently imply different observa- 
tional properties. 

A well-known property of the climate sys- 
tem is its chaotic nature, arising from the com- 
bined and concomitant action of different non- 
linear feedbacks, that may give rise to hystere- 
sis cycles and multi-stable solutions. This im- 
plies that, as occurred several times to Earth, 
a planet could be habitable and inhabited even 
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Fig. 1. Habitability map on the semi-major 
orbital axis vs eccentricity plane (color-code 
in the lower bar; the habitability criterion 
is the "complex" habitability as defined in 
Silva et al. (2017b). Those solutions admitting 
climate warm-snowball bistability are shown 
with filled circles. Their color-code in the 
right bar gives the snowball probability. Most 
bistable and habitable solutions are found 
within the same parameter range. See Murante 
et al. (2020) for more details. 


if in a snowball state. This possibility can be 
accounted for with a probabilistic approach, as 
we have shown in Murante et al. (2020, M20 
hereafter). 

By analyzing a few x10^ simulationg?]|M20 
explored the frequency and parameter depen- 
dence of bistable warm-snowball climate so- 
lutions. The considered parameters — pressure, 
semi-major axis, eccentricity, obliquity — were 
coupled to a range of different initial condi- 
tions (represented by the initial temperature) 
to run the simulations. M20 found a probabil- 
ity x 10% of bistable solutions within the ex- 
plored range of parameters, but interestingly 
found an indication of a connection between 
habitability and bi-stability, i.e., that the range 
of parameter values giving rise to bistability 
is similar to range allowing habitable solutions 


(Fig.[I). 


? The simulations database, ARTECS, is avail- 
able at https://wwwuser.oats.inaf.it/exobio/climates/ 
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Fig. 2. Global mean surface temperature of the 
Earth during the Archean. The colored squares 
indicate the results obtained in the present 
work for the atmospheric compositions in the 
legend (med/low CO» is for pco,=0.04/0.01 
bar; high/med CH, is for 1000/100 ppm). 
The gray curves are the lower and upper lim- 
its of the temperature constraints obtained by 
Catling & Zahnle (2020), while the magenta 
upper limit and the cyan circle represent ge- 
ological constraints on temperature by Hren 
et al. (2009) and Blake et al. (2010), respec- 
tively. The dashed horizontal line is for the wa- 
ter freezing point. 


In addition to the ice-albedo feedback, we 
have recently introduced in ESTM also the 
vegetation-albedo feedback (Bisesi et al., in 
prep.). In fact, the temperature dependence of 
the vegetation coverage of lands heavily affects 
the albedo, being significantly smaller for veg- 
etation rather than for bare continents (~ 0.1 
vs ^ 0.3) On Earth, the vegetation/albedo 
feedback is known as Charney mechanism 
(Charney et al. 1975). It affects the Earth cli- 
mate and, depending on the planetary proper- 
ties, it can increase the width of the circum- 
stellar HZ by up to 10% on its outer border. 
The Charney mechanism may also increase the 
global habitability of the planet by a similar 
amount. 


3.2. The Archean Earth 


The Archean eon spans the time from 4 to 2.5 
Gyr ago, i.e. roughly 1/3 of the Earth evolu- 
tion, and it probably accommodates the emer- 
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gence of life. It is therefore particularly im- 
portant to try to understand the evolution of 
the environmental conditions ruling habitabil- 
ity during such a long period of time. This eon 
was characterized by an anoxic atmosphere, a 
shorter rotation period, a very low coverage 
of continents, and by a fainter insolation (~ 
0.7 — 0.8 the current value, due to stellar evolu- 
tion). In fact, the well-known Faint Young Sun 
problem refers to trying to reconcile the geo- 
logical evidences of the presence of liquid wa- 
ter during the Archean with such a low energy 
input from the Sun. Among several co-factors, 
such evidence is generally ascribed to larger 
amounts of mainly CO, and/or CH4, able to 
provide a greenhouse forcing much larger that 
the current value (i.e. > 30 °C). According to 
several geological constraints (e.g. Catling & 
Zahnle 2020) the Earth planetary/atmospheric 
properties have evolved during the eon, in such 
a way to maintain on average a temperate sur- 
face. 

EOS-ESTM is particularly suited to study 
the large parameter space that may have al- 
lowed habitability of the Archean Earth. We 
present below some preliminary results that 
we have performed by considering a range of 
atmospheric compositions (pco,—0.01 to 0.1 
bar, and CH4 abundance = 0 to 1000 ppm, 
for a total N5-dominated dry pressure of 1 
bar with 60% relative humidity), and at the 
same time by accounting for the time evolution 
of the luminosity of the Sun (Gough 1981), 
of the rotation period (Bartlett & Stevenson 
2016), and of the fraction of land coverage 
(Collerson & Kamber 1999). For more details 
see Bevilacqua (2022); Simonetti (2022). 

In Fig. |2} we show the evolution of the 
global mean surface temperature estimated 
from available data and compared to our sim- 
ulations. The gray curves are the lower and 
upper limits of temperatures with 9596 confi- 
dence intervals obtained by Catling & Zahnle 
(2020). The magenta upper limit and the cyan 
circle are the constraints from geological data 
provided respectively by Hren et al. (2009) and 
Blake et al. (2010). The squares with error 
bars are the results obtained from our simula- 
tions. These are computed as follows. We have 
split the 1.5 Gyr-length of the eon into three 


500 Myr-long intervals. For each step, we se- 
lected among our currently computed simula- 
tions (that at each age adopt the appropriate 
insolation and planetary properties) a compo- 
sition able to provide a habitable surface tem- 
perature. Specifically: (a) blue square, 4-3.5 
Gyr, pco,-0.04 bar, CH4-1000 ppm; (b) red 
square, 3.5-3 Gyr, pco,=0.04 bar, CH4-100 
ppm; (c) green square, 3-2.5 Gyr, pco,=0.01 
bar, CH4-100 ppm. The temperatures we ob- 
tained are consistent with the geological con- 
straints derived by Hren et al. (2009) and Blake 
et al. (2010). 

In each phase, by considering the Earth 
as a benchmark for exoplanets, the observable 
spectral features would be different. We are 
working to produce a finer model grid in or- 
der to explore a larger range of possible so- 
lutions, and also to simulate for each case, 
the corresponding expected observable atmo- 
spheric spectrum. 


4. Connecting models of climate and 
habitability to observations 


The atmospheric characterization of rocky 
exoplanets may be hopefully achievable al- 
ready during the JWST lifetime for M-type 
hosts, but it is considered a fundamental tar- 
get for already planned or under design space- 
and ground-based instruments during the next 
decades (e.g. ARIEL, ELT, LUVOIR, HabEX, 
LIFE). Depending on the line-of-sights con- 
figurations, orbital distances and wavelengths, 
spectra will be collected via transmission, re- 
flection or emission. The same atmospheric 
structure and composition that dictates the cli- 
mate and the surface environmental conditions 
for habitabilityP] rules also the spectral emis- 
sion that will be measured. It is therefore nec- 
essary to model at the same time climate and 
habitability, with the associated atmospheric 
spectrum. 

As an example of the connection between 
observables and habitability, in Fig. B]we show 


? Habitability can be quantified with different 
possible temperature-based criteria once the climate 
state has been computed. Also, it can be extended to 
include the protective role of the atmosphere against 
ionizing radiation (see Silva et al. 2017b). 
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Fig. 3. Habitability maps in the atmospheric mass (log p/g in gr/cm?) vs insolation plane, as in 
Silva et al. (2017b), with superimposed the contour lines of constant OLR (left panel) and TOA 
albedo (central panel). These curves decouple and intesect at large enough p/g (right panel, 
where for clarity only one of the albedo lines is shown in white together the OLR lines). Since 
OLR and TOA albedo, together with insolation, are observables, they could provide a measure 
of the habitability. The red-yellow circle is the position of the Earth. The thick red line marks the 
atmospheric mass limit below which the surface radiation dose of secondary particles of Galactic 
Cosmic Rays would be > 100 mSv/yr (see Silva et al. 2017b). 


the habitability map as a function of atmo- 
spheric mass (p/g) and insolation, as in Silva 
et al. (2017b), with superimposed the cor- 
responding contours of constant OLR (left 
panel), and TOA albedo (central panel). At 
large enough p/g values, these contour lines 
decouple and intersect. The observational eval- 
uation of OLR and Bond albedo, together with 
the measured insolation, could therefore in 
principle provide an estimate of the habitabil- 
ity. This sort of statistical inferences requires 
large numbers of model simulations, which we 
can perform with our modelling tools. 


4.1. Transmission spectra: effect of 
refraction 


Exoplanets discovered with the transit method 
— so far the majority — are in principle good 
candidates for measuring their atmospheric 
spectrum in transmission. However, transmis- 
sion spectra are affected by refraction, which 
impacts the accessible atmospheric layers. The 
deflection angles depend on atmospheric com- 
position and density, and on the angular size of 
the host star as seen from the planet. Therefore 
refraction is an issue for the relatively dense at- 
mospheres of terrestrial-type rather than giant 
planets, and in particular for habitable rocky 


planets around G-type stars due to the implied 
geometry. 

We present preliminary results from Maris 
et al. (in prep.) in simulating observations 
that account for refraction. We highlight that 
the same atmospheric models that, with EOS- 
ESTM, determine the climate state, are then 
used to compute the associated refraction and 
observability of the spectral features. 

Fig. |4| left panel, shows the fraction, with 
respect to no refraction, of the integrated band 
extinction for 3 representative bands for O5 
(0.759 — 0.775 um), CO» (1.565 — 1.591 um) 
and H5O (0.860 — 1.020 um), that would be 
Observed, for a model of the modern Earth, at 
different sampled atmospheric heights. For this 
model, only the layers above 11.8 km (i.e. in 
practice above almost all the troposhere con- 
taining most atmospheric mass) would be al- 
ways visible during transit, implying that the 
depths of the considered bands would be re- 
duced to ~ 60% for O and 40% for CO2 
compared to their full depth. Instead the wa- 
ter vapor band would be totally missed. Note 
that the shown example is for the best ideal 
clear-sky case. In the right panel we show an 
illustrative simplified example of the possible 
effects of clouds on the transmittance of the 
molecular absorption spectra. The blue and or- 
ange spectra are without refraction, but for the 
latter spectrum a 10 km cut has been applied 
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Fig. 4. Left: Fraction of observed integrated band extinction with respect to no refraction, vs the 
minimum probed atmospheric height, for three representative bands for O, (0.759 — 0.775 um), 
CO, (1.565 — 1.591 um), H20 (0.860 — 1.020 um). This example is for a current Earth-like 
atmosphere in the best case of clear-sky conditions. The vertical dashed and solid lines show the 
minimum heights probed at transit start and center, respectively (7.4 and 11.8 km for this model). 
All layers above the latter minimum height are always visible during transit. Right: transmittance 
spectrum of molecular absorption in the ~ 1 — 2 um region. The blue and orange lines are for no 
refraction, but in the latter case a 10 km cut is applied to represent a 10096 cloud coverage. This 
case is almost superimposed to the green line, that is the transmittance that is obtained including 
refraction. See text for more details. From Maris et al. (in prep). 


in the sampled atmospheric layers, in order to 
schematically represent a 10096 cloud cover- 
age. This case is almost superimposed to the 
green line, the transmittance that is obtained 
accounting for refraction. Note that the shown 
examples do not include the Rayleigh scatter- 
ing. This and other effects and complications 
will be presented in Maris et al. (in prep.). 


While clouds (and hazes) will be a se- 
vere issue for the detection of spectral features 
and biosignatures also for habitable planets in 
close orbits around M-stars hosts, refraction 
will severely limit the possibility of detecting 
features precisely for habitable terrestrial ana- 
logues in any case. In these situations, emis- 
sion spectra could provide better probes. 


4.2. Emission spectra 


Emission/reflection spectra are valuable tools 
that both allow for the characterization of the 
disk-integrated vertical pressure-temperature 
structure of the atmosphere and give informa- 
tion on the horizontal heat distribution. An 
example of this type of analysis has been 


conducted on WASP-103b (Kreidberg et al. 
2018). Currently, a few low resolution emis- 
sion/reflection spectra of hot and ultra-hot 
Jupiters are available, while proposed future 
missions such as the Large Interferometer For 
Exoplanets (LIFE, Quanz et al. 2022) will fo- 
cus on the study of rocky, temperate planets. 
As an intermediate step in producing the 
OLR and TOA albedo lookup tables, EOS cal- 
culates the synthetic emission/reflection spec- 
trum of the planet. An example of this is shown 
in Fig. |5| which has been produced using an 
Earth-like atmosphere containing N2, O2, 360 
ppm of CO», 1.8 ppm of CH, and a varying 
amount of H5O, calculated fixing the relative 
humidity to 60%. The surface temperature has 
been fixed to 280 K and the vertical pressure- 
temperature profile has been specified in input, 
using a moist adiabatic lapse rate for the tro- 
posphere and a 200 K isothermal stratosphere. 
The surface albedo has been set to 0 and the in- 
coming stellar radiation has been produced by 
a blackbody at 5780 K, with a total flux nor- 
malized to the Earth’s Solar Constant (1361 W 
m ?). The calculation has been carried on in 
the [0.34, 200] um range at a A/A4 resolution 
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Fig. 5. Top of atmosphere SED for an Earth- 
like planet with a surface temperature of 280 K, 
a tropopause at 200 K and a surface reflectance 
of 0. The instellation is produced by a black- 
body at 5780 K, emulating the Sun. The red 
line represents the OLR only, while the yel- 
low line represents the sum of OLR and the 
reflected stellar radiation due to atmospheric 
Rayleigh scattering. The dotted lines describes 
the blackbody curves for a 280 K and a 200 K 
emitter. 


of 3000. This example, calculated for a sin- 
gle atmospheric column, shows that, with the 
disc-integrated emission spectrum appropriate 
to each EOS-ESTM simulation, we will be able 
to test the detectability of habitable exoplanets 
with future space missions (e.g. LIFE). 
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